2D Metal-Insulator transition as a percolation transition 
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By carefully analyzing the low temperature density dependence of 2D conductivity in undoped 
high mobility n-GaAs heterostructures, we conclude that the 2D metal-insulator transition in this 
system is a density inhomogeneity driven percolation transition due to the breakdown of screening 
in the random charged impurity disorder background. In particular, our measured conductivity 
exponent of ~ 1.4 approaches the 2D percolation exponent value of 4/3 at low temperatures and 
our experimental data are inconsistent with there being a zero-temperature quantum critical point 
in our system. 
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Ever since the pioneering observation by Kravchenko 
et al. 1] of an apparent two-dimensional (2D) metal- 
insulator transition (MIT) in Si-MOSFET inversion lay- 
ers, the nature of the transition has remained a contro- 
versial enigma. The 2D MIT has been claimed by some, 
most notably by Kravchenko and collaborators 1,^, to 
be an interaction driven and carrier density tuned T=0 
quantum phase transition whereas others have argued 
|3 that it is a disorder-driven crossover phenomenon be- 
tween weakly and strongly localized 2D electron states. 
Density dependent scanning studies of the 2D chemi- 
cal potential |j|, surface acoustic wave measurements on 
low density 2D electrons |3, and compressibility measure- 
ments la, I3, as analyzed by direct numerical simula- 
tions of the energetics of the 2D disordered system, 
indicate very strong density inhomogeneities in the 2D 
system around the critical density (ric) for the putative 
2D MIT. Also, the critical density Uc is highly system- 
dependent and varies very strongly with the impurity 
disorder in the 2D system. In relatively highly disor- 
dered Si MOSFETs the typical n^ « lO^cm-^ and in 
very high quality 2D n-GaAs heterostructures (the sub- 
ject matter of our work presented in this paper) the typ- 
ical Uc ~ lO^cm^^. The strong dependence of Uc on 
disorder, the long range nature of the charged impurity 
disorder potential in 2D semiconductor structures, the 
observed density inhomogeneities around ric, all taken 
together, suggest a percolation-type transition underly- 
ing the 2D MIT. Indeed such a percolation transition for 
2D systems was discussed by Efros fifteen years ago [2, 
and has recently been suggested in the context of the 2D 
MIT phenomena by several authors |1 ITollllllTllT3l| . 
In this Letter we provide rather compelling experimental 
evidence in support of the 2D MIT being a percolation 
transition, at least for our experimental system. 

The basic picture of the percolation transition for 2D 
MIT is simple and highly physically motivated. As the 
carrier density (n) is lowered in a 2D system, screen- 
ing becomes progressively weaker and strongly nonlin- 
ear. A small decrease in n leads to a large decrease in 



screening, and eventually to a highly inhomogeneous 2D 
system as the electron gas at low enough carrier den- 
sities becomes unable to screen the disorder potential 
with the individual charged impurity scattering poten- 
tials getting unscreened and their random distribution in 
the 2D system. This gives rise to the well-known ran- 
dom "hill-and-valley" potential landscape with the 2D 
carriers repelled from the potential hills and accumulat- 
ing at the potential valleys in contrast to the high carrier 
density homogeneous situation. Once these "depletion 
(or denuded) zones", associated with the disorder po- 
tential hills are numerous enough to prevent percolating 
conducting paths to span the 2D system connecting the 
"valley" regions, an effective 2D MIT transition takes 
place with the system being an effective metal (insula- 
tor) for n > (<) Tic where Uc is the critical percolation 
density. This percolation picture is particularly germane 
to 2D semiconductor systems because the disorder po- 
tential here arises from the presence of random charged 
impurities in the system making electronic screening a 
key ingredient in the effective disorder seen by the car- 
rier system. Such a density-driven percolation transition 
is the hallmark of the long-range Coulomb disorder po- 
tential in the system since the Coulomb disorder can be 
tuned by changing the carrier density through the non- 
linear screening mechanism. The key qualitative features 
of the 2D MIT phenomenon, namely the correlation be- 
tween Tic and the sample quality as well as the density 
inhomogeneity around Uc, are entirely consistent with a 
percolation transition in the screened Coulomb disorder 
potential. 

Motivated by these percolation considerations we have 
carefully analyzed the density dependent conductivity 
a{n) in ultra- high mobility n-GaAs undoped heterostruc- 
tures which have demonstrated [14] the lowest value of ric 
(- 2 X lO^cm-2) foj. 2D MIT so far in the literature. Most 
of the 2D MIT literature has emphasized the strong tem- 
perature dependence in the effective metallic (n > ric) 
phase, which is now well-understood [13 as a manifesta- 
tion of the strong temperature dependence of the effective 



Coulombic disorder seen by the low-density 2D carriers 
through the temperature dependence of screening. In 
this Letter, we analyze the power law behavior of cr(n) 
in the context of density dependent screening and a per- 
colation transition. We have also carried out a detailed 
Boltzmann theory calculation |1?t| using realistic quasi- 
2D parameters appropriate for our system, assuming the 
carrier transport to be dominated by charged impurities 
distributed randomly at the interface - the main physics 
being the temperature and wave vector dependent linear 
screening of the long-range Coulomb scattering potential 
of the charged impurities. We have earlier shown [lj| that 
such a screening description provides an excellent qual- 
itative agreement with experiment for the temperature 
dependent transport properties of the effective metallic 
phase in our system. 

The two samples used in this study are undoped 
GaAs/AlGaAs heterojunctions with field-induced carri- 
ers. From the surface, the epilayers consist of 60 nm n-f 
GaAs followed by a 600 nm Alo.3Gao.7As barrier, and 
then 1000 nm of GaAs. The 2D electrons are induced 
from the ohmic contacts with a density that is propor- 
tional to the voltage on the n+ GaAs layer that serves as 
a top gate. Standard four terminal lock-in measurements 
are made to measure the resistance over a wide range of 
density. Conductivity of these square samples is deter- 
mined using van der Pauw techniques. In Sample A the 
density spans a range of 0.15 to 7.5 x lO^'^cm"^ with a 
very high maximum mobility of 8.5 x lO^cm^/Vs at high 
density and sample B spans a density range of 0.5 to 
13 X lO^^cm-^ with a lower mobility of 3.3 x lO^cm^/Vs 
at high density. 

The experimentally measured conductivity cr(n) is 
shown in Fig. 1 along with the calculated ^.13| fj{n) curves. 
We first discuss a few salient qualitative features of the 
density dependent transport properties of our results. At 
high densities {n > lO^'^cm^^), the conductivity depends 
on the density approximately as cr '^ n" with a w 1.6. 
We emphasize that this is not a strict power law since 
a = a{n) depends weakly on the density. This high 
density behavior is completely consistent with our Boltz- 
mann theory based calculations assuming the conductiv- 
ity being limited by linearly screened charged impurity 
scattering. As n decreases, a starts decreasing faster with 
decreasing density and the experimental conductivity ex- 
ponent a(n) becomes strongly density dependent with 
its value increasing substantially above the high-density 
value of a w 1.6. At the lowest density for sample A, 
a « 5. Although a part of this strong density dependence 
of a{n) at low density can be understood as arising from 
the strong suppression in screening at low carrier densi- 
ties, we find that screening in a homogeneous electron 
gas fails qualitatively in explaining the a{n) behavior 
at low densities whereas it gives quantitatively accurate 
results at high densities. As has been found from di- 
rect numerical simulations % ;9| , homogeneous screening 
of charged impurity disorder breaks down at low carrier 
densities with the 2D system developing strong inhomo- 
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FIG. 1: Experimentally measured (symbols) and Boltzman 
scattering calculated (lines) values of conductivity for sample 
A (la) and sample B (lb). 



geneities leading to a percolation transition at n = ric- 
For n < Tic, the system is an insulator containing isolated 
puddles of electrons with no "metallic" conducting path 
spanning through the whole system. 

The percolation scenario also naturally explains the 
qualitative aspects of the observed strong temperature 
dependence of the effective metallic phase in low-disorder 
high-quality samples for n > ric in the 2D MIT phe- 
nomenon. In particular, low disorder ensures a low 
value of Tic, which then automatically leads to a strongly 
temperature dependent screening correction producing a 
temperature dependent conductivity as both T /Tp and 
qxp/'^kp are effectively large for low densities[l5l|. In 
samples with large values of n^ (low-quality samples) the 
temperature dependent screening effects are weak in the 
effective metallic phase [n > ric) since the carrier den- 
sities are high with the associated qxp/'^kp and T/Tp 
being "small" [13. 

In Fig. 2 we fit our measured conductivity a{n) to the 
expected [9llla| percolation 'critical' behavior 



a{n) = A{n - Uc)^ 



(1) 



where 6 is the conductivity percolation exponent char- 
acterizing the vanishing of the conductivity. Nonlinear 
curve fitting to Eq. (1) with proper weighting of the 
conductivity is performed for a range of the low density 
data to determine values of S, ric and A. For sample A 
the data used in the fit is 1 x 10^ < n < 4 x lO^cm"^, al- 
though the results do not depend sensitively on the exact 
range of data used in the fit. Note that for T = 0.05 K 
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FIG. 2: Experimental conductivity (symbols) fit to Eq. 1 
(lines) for both Sample A and B as indicated in the legend. 
The same results are displayed as a function of n (2a) and 
n — He (2b) to better show the low density region. 



and n < 2.3 x lO^cm"^ a large lock-in quadrature signal 
prevents measurement of the conductivity. For Sample 
B, the 'goodness' of the fit is more strongly dependent on 
the range of data used for fitting. For the data shown, 
the lowest density point is excluded from the fit, and 
the fitting range is 5 x 10^ < n < 8 x lO^cm^^. The 
values of the relevant parameters are 5 = 1.4 ± 0.1 and 
Uc = 0.18 ±0.01 X 10 



^°cm ^ at T = 50 mK for sample 



A and 1.5 ± 0.1 and n^ = 0.28 ± 0.02 x lO^^cm-^ for 
sample B. These values of S are close to the known 2D 
percolation exponent of 4/3,16] ■ 

In Fig. 3 we show the dependence of the critical ex- 
ponent S{T) and the critical density nc{T) of sample A. 
At low temperature, the experimental conductivity ex- 
ponent 6(T) approaches a value S « 1.4 (1.5 for sample 
B), which is close to the expected 2D percolation expo- 
nent. The decrease of S with T for T > 0.2K is caused 
by the increasing importance of phonon scattering which 
should be disregarded for our discussion of the 2D MIT 
itself. The increase of 6 (for T < 0.2 K) and decrease of 
ric with increasing temperature can be well-understood 
by adding a non-critical temperature dependent contri- 
bution /(T), where /(T ^ 0) = and /(T) ^ for 
T 7^ 0, to the right-hand side of Eq. 1 and realizing 
that at any finite T the system, by definition, is an "ef- 
fective" metal since the conductivity is non-zero at all 
finite temperatures. This automatically implies that the 
effective 5{T) increases with T and the effective nc{T) 
decreases with T, since 5 and ric are strictly defined only 
at T = and at all finite T the system tends to behave 
as an effective metallic phase - in fact, for sufficiently 
high values of T (> li^) one cannot distinguish at all 
between the metallic and the insulating phase by look- 
ing at the conductivity. The fact that critical behavior 
manifests over about a decade of density and two decades 
of conductivity as well as our finding of a conductivity 
exponent consistent with the 2D percolation transition 
suggests that the 2D MIT in these low disorder n-GaAs 
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FIG. 3: Parameters 5 and Uc for sample A for a range of 
temperature. Error bars are the 95% confidence interval as 
determined from the nonlinear fitting shown in Fig. 2. Sample 
B shows much weaker temperature dependence of 5 and Uc in 
this range. 



systems is a percolation transition. 

It is important to point out that the critical density 
ric we obtain by the optimal fitting of experimentally 
measured a{n) to the conductivity scaling formula of Eq. 
1 is consistently lower than the corresponding estimate 
{n'c) based on the sign of da/dT as is often done in the 
literature to distinguish an effective metal from the in- 
sulator. This is perfectly understandable based on the 
"qu antum-classical crossover" scenario discussed in ref. 
llSl - at low carrier densities the effective metallic phase 
exhibits an insulating da/dT (> 0) sign down to rather 
low temperature, and therefore any determination of the 
critical density n'^ based purely on the sign of da/dT 
will produce n'^ > ric- Thus, the 2D MIT transition 
density cannot be ascertained by looking for the point 
da/dT\n=n' = - the sign of da/dT is not a good indi- 
cator for the critical density. Our analyses as shown in 
Figs. 1-3 imply that the transition itself is better under- 
stood by concentrating on the density dependence oia{n) 
at fixed (and very low) values of temperature whence the 
2D MIT can be studied approaching entirely from the 
metallic regime with the sign of da/dT playing no role in 
the analysis. The effect of electron-electron interaction 
(beyond just the nonlinear screening mechanism leading 
to the electron 'puddles' or 'droplets') is mainly to de- 
crease the effective value of the critical density at which 
the percolation transition occurs, since the main effect 
of interaction is to 'homogenize' the density of the elec- 
tron liquid in the process reducing the effect of density 
inhomogeneities arising from the random charged impu- 



rity centers. Without interaction effects, percolation will 
occur at higher carrier densities Q- 

Before concluding we emphasize that one obvious and 
immediate consequence of the 2D MIT being a disorder 
driven percolation transition in our system is that it is 
manifestly not a quantum critical phenomenon as has 
sometimes been suggested in the literature. We have, in 
fact, attempted to fit our density and temperature depen- 
dent conductivity cr(ri,T) to the usual quantum critical 
scaling form P,[|, a{n,T) - F{T/To), where F is the 
universal scaling function with Tq = To(n) ~ jn — nd''' 
where the exponent j = vz is a product of the corre- 
lation exponent v and the dynamical exponent z. We 
have not been able to obtain any kind of scaling col- 
lapse of our a{n,T) data as a function of T/Tq. The 
inconsistency of our results with quantum criticality is 
obvious from the percolation scaling fits to Eq. (1) that 
our data exhibit so well. The fact that our fitted per- 
colation exponent S and critical density ric both depend 
on temperature already indicates an absence of quan- 
tum critical scaling in the data. A more subtle point 
is that the usual quantum critical 'fan' diagramjl^ indi- 
cates that the scaling regime should expand with increas- 
ing temperature (at low enough temperatures) since the 
quantum critical regime spreads out in density (around 
Uc) at higher values of T. The expected quantum critical 
'fan' behavior is in strong disagreement with our observed 
behavior of a{n,T) where we find that the best scaling 
of (T(n) ~ (n — ficY occurs at our lowest measurement 
temperatures, and the scaling regime (i.e. the regime in 
carrier density) definitively and systematically shrinks as 
temperature increases in sharp contrast to the quantum 
scaling fan behavior. It cannot, of course, be ruled out 
that the quantum critical regime ends at temperatures 
way below our lowest measurement temperature (^ 50 
mK). We believe that our data and the percolation anal- 
yses presented in this paper rules out a quantum critical 
phenomenon underlying 2D MIT. 

We conclude by emphasizing that we provide in this 
Letter compelling evidence involving both experimental 
transport data and detailed theoretical analyses in sup- 
port of the 2D MIT phenomenon being a 2D percolation 
transition in our system where the carrier system breaks 
up into an inhomogeneous collection of disconnected local 
'droplets' of electrons localized in the potential valleys of 
the long-range Coulombic disorder potential. Our mea- 
sured conductivity exponent 5 w lA — 1.5, approaches 
the expected 2D percolation exponent of 4/3, and we see 
impressive scaling behavior in our data consistent with a 
percolation transition. The percolation transition occurs 
at a critical density ric considerably lower than the corre- 



sponding 'critical' density nj, one would infer for 2D MIT 
based just on the temperature dependence of ct(ti, T). By 
comparing with our theoretical calculations, we establish 
that the effective disorder underlying 2D MIT is the scat- 
tering from screened charged impurity centers randomly 
distributed near the semiconductor-insulator interface. 
As screening weakens at low carrier densities, the sys- 
tem breaks up into inhomogeneous carrier droplets sep- 
arated by denuded or depleted zones, leading eventually 
to the percolation transition as was discussed by Efros 
some years ago H • Our transport data validates the ex- 
istence of a percolation transition in 2D MIT but cannot 
really provide a qualitative description of the underlying 
percolation process - for example, we cannot rule out per- 
colation scenarios which are somewhat different from the 
physically appealing (but perhaps somewhat simplistic) 
picture of Efros where the percolation transition occurs 
primarily through the depleted zones becoming numer- 
ous enough at low enough carrier densities. More ex- 
perimental work, involving simultaneous measurements 
of compressibility, capacitance, transport, local chemical 
potential, and perhaps even local conductivity, is clearly 
required to better understand the details of the percola- 
tion mechanism and the nature of the inhomogeneous 2D 
carrier system. Our work also brings up important ques- 
tions regarding the universal nature of 2D MIT itself - for 
example, we urge experimental analyses of <y{n) at fixed 
low temperatures in other 2D systems to check whether 
this percolation scenario is universal or just applies to our 
ultra-high mobility 2D n-GaAs system. In this context 
it is perhaps quite interesting to point out that even the 
extensively studied 3D MIT in doped Si systems has been 
speculated [ig to be a possible inhomogeneity driven per- 
colation transition (similar to what we find for 2D MIT 
in our work) , and such a percolation transition may very 
well be a generic mechanism for metal-insulator transi- 
tions in semiconductors where screened charged impurity 
potential is the dominant disorder mechanism. 
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